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It is at present well known that conditions at some massive binary systems allow accel- 
eration of particles and production of the GeV-TeV 7-rays. However, which particles are 
responsible for this emission and what radiation processes are engaged is at present not 
completely clear. We discuss what parameters can determine the acceleration process of 
particles and high energy radiation produced by them within massive binary systems. 
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1. Introduction 

The high energy 7-ray emission from the binary systems containing energetic pulsars 
has been considered since late 70-ties (Cyg X-3L^or LS I +61 303 El). However, the 
first enthusiastic reports on the positive 7-ray signals in the TeV-PeV energies from a 
few binaries (e.g. Her X-l, Vela X-l, Cyg X-3, ...) have not been confirmed by more 
reliable observations with the Whipple Cherenkov telescope SI. These likely false 
reports initiated theoretical investigation of possible propagation of 7-rays injected 
into the volume of the binary system s it; Ua tj 0I1 changed with the launch of 

the EGRET telescope on board of Compton GRO. A few binaries have been found 
in the large error boxes of the EGRET GeV 7-ray sources (e.g. LS 5039 ^, Cyg 
X-3 13, LS I +61 303 9 , Cen X-3 10 ). These claims initiated more detailed studies 
of the production and propagation of 7-rays within the binary systems such as, 
e.g. simulations of the 7-ray pr oductio n in the anisotropic Inverse Compton e pair 
cascade in the stellar radiation I H I -^ I 

The breakthough came with the discoveries of 7-ra y bi naries at TeV energies 
by the modern Cherenkov telescopes (LS2883/PSR1259Q31, LS 503912, L $ I +61 
303 These discovaries initiated efficient investigation of different scenarios for 
the 7-ray production in the massive binary systems. In this paper I discuss what 
parameters of the binary systems can influence the acceleration of particles and 
subsequent production of GeV- TeV 7-rays in these systems. 
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2. Gamma-ray observations of binary systems 

The first binary system containing an energetic pulsar LS2883/PSR1259 was de- 
tected at the TeV energies to be related to the crossing the periastron passage ^^1. 
The TeV emission seemed to be related to the moments of the intersection of the 
pulsar through the equatorial wind of the massive star (B e ty pe) in this system. 
GeV emission has been also discovered from this system showing an unex- 
pected brightening at about a month after the periastron passage On the other 
hand, LS 5039 and LS I +61 303 (also now suspected to contain pulsars due to the 
characteristic radio morphology, Dhawan et al. and Ribo et al. showed clear 
modulation of the TeV signals with the periods of the binaries with the ma ximum 
emission when the compact objects were not far from the inferior conjunction EOEl] 
A clear correlation of the TeV and X-ray emission has been reported in the case 
of LS I +61 303 an d LS 5 039, thus supporting their production by this same pop- 
ulation of electrons l ^ l ^ l Q e y -y-ray emission from t hese two binaries show 
clear unticorrelation with the TeV 7-ray light curve l ^4 | 25[ j s k nown a t present 
that the TeV emission from at least from LS I +61 303 shows long s cale variability 
during which the 7-ray light curve changes significantly its shape idBlALi The GeV 
and TeV spectra of LS 5039 and LS I +61 303 show two di stinct components with a 
transition in the energy range between 10-100 GeV^ES] Recentl y, two additional 
massive binary systems, HESS J0632+057 detected in TeV energies l 28 l 29 l an d 1FGL 
J1018. 6-5856 detected in GeV energies^!, show modulation of the 7-ray signal with 
the periods of the associated binary systems. 

Another very compact binary system, CygX-3, belonging to the microquasar 
class, has been also detected in GeV 7-rays 1 1 I The 7-ray emission observed from 
this system is correlated with the major radio flares ^5 The emission shows a modu- 
lation with the period of the binary system, with the maximum close to the location 
of the compact object behind the companion star. Up to now, TeV 7-ray emission 
has not been discovered from Cyg X-3, in spite of extensive observations with the 
MAGIC telescope^S. In the case of another microquasar, the massi ye bina ry system 
Cyg X-l, there are evidences of transient GeV- TeV 7-ray emission I 34 | 35 | T he flare 
of TeV 7-ray emission, almost instantaneous with a flare in the hard X-rays I 34 | 36 |^ 
has been reported from Cyg X-l. It lasted for less than an hour, around the su- 
perior conjunction of the compact object. Also a separate flare of 7-ray emission 
above 100 MeV, on a time scale of 1-2 days, has been reported recently by the Ag- 
ile Team EH. These GeV 7-ray flares seem to be rather exceptional during the 300 
day observation period with the Agile telescope The 7-ray emission from other 
microquasars within binary systems has not b een di scovere d up to n ow in spit e of 
extensive observations, see e.g. GRS 1915+105 EHEH^ SS433 | 40 | ^ 1 | 2 Sor Sco X-l^. 

Only one massive binary system containing two massive stars, Eta Carinae, 
has been l ikely dete cted by the Agile and Fermi-LAT telescopes at energies below 
^100 GeVMESl. The high energy emission of Eta Carinae is characterised by an 
intriguing two-component spectrum which can be approximated by two power laws, 
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Fig. 1. The geometry of the interaction of relativistic particles with the radiation field within the 
binary system. Particles are accelerated at shocks within the jet or at the shock wave which results 
from the pulsar wind-stellar wind collision. The geometry of the acceleration and the interaction 
process may/or may not be similar, depending on the parameters of the binary system. 

first one with an exponential cut-off at a few GeV and the second extending up to 
^100 GeV. It has been reported that the second 7-ray component varies with the 
phase of the binary system, thus supporting the evidence that it has to come from 
within the binary system ^1 It will be very interesting to search for the possible 
TeV emission from this binary system. 

3. Conditions within massive binary systems 

7-ray emission has been detected from two types of massive binary systems. In 
the first one it is expected that the matter accreting onto a black hole creates 
an accretion disk. From the inner part of the disk, a jet is launched similarly to 
the process observed in active galactic nuclei (the so-called microquasar scenario). 
Particles are expected to be accelerated in the jet, which is usually considered to 
propagate perpendicularly to the disk plane and to the plane of the binary system 
(see Fig. 1). In the second case, the binary system contains either two massive 
stars or an energetic pulsar and a massive star. Both stars produce strong stellar or 
pulsar winds. Particles are expected to be accelerated at the shock structure which 
appears as a result of the collisions of the winds. They can be also injected from 
the inner magnetosphere of the pulsar. The geometry of injection of relativistic 
particles into the binary system in these two scenarios may (or may) not differ 
significantly. If the pressure of both winds in the binary system is comparable, then 
the shock structure roughly has a plane shape. However, in the case of the winds 
with different strengths, the shock can likely bend around one of the stars producing 
different injection geometry for relativistic particles. 

The radiation process by accelerated particles within the binary system can also 
occur in a complicated environment. After all, the magnetic field around the massive 
stars can have a complicated structure (see Fig. 2). Depending on the strength of the 
wind and the stellar rotation velocity, the magnetic field structure is well described 
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Fig. 2. The magnetic field around massive stars have a very complicated structure (left figure), 
close to the star it is dipolar, at larger distances, it is radial, and far away from the star, it 
becomes toroidal. The winds of Be type stars are strongly inhomogeneous, being slow and dense 
in the equatorial plane of the rotating stars and fast and rare in the poloidal region (right figure). 



by the dipole approximation (B(r) oc r~ 3 ). This usually happens only in the region 

very close to the stellar surface (up to 2 stellar radii in the extreme case). At 

farther distances, the magnetic field becomes radial (B(r) oc r~ 2 ), and finally far 

away from the star, it is expected a toroidal structure (B(r) oc r _1 ). The transition 

between these two last regions is typically at the distance of ~ 10 stellar radii. It is 

l47l 

determined by the rotational velocity of the star l - iiJ . 

Also the winds produced by the massive stars are not uniform. For example, 
the Be type stars (which are expected to be present in LS2883 and LS I +61 303) 
have dense and slow equatorial winds and fast and rare poloidal winds (see Fig. 2). 
Therefore, it is expected that the shock structure within the binary system and the 
conditions for acceleration of particles can change (sometimes even very abruptly) 
with the orbital phase of the binary system. 

Massive stars in 7-ray compact binary systems produce a very strong radiation 
field. 7-rays with sufficiently large energies are efficiently absorbed on this radiation 
field. In fact, detailed calculations of the optical depths for 7-rays are clearly above 
unity for the parameters typical for the massive stars These optical depths 
strongly depend on the location of the injection place within the binary system 
and on the propagation angles in respect to the direction towards the massive star 
(note that for very compact binary systems the dimensions of the stars have to be 
taken into account). Due to the similarities of the cross sections for the 7 — 7 — > e 
absorption process and the Inverse Compton (IC) scattering of soft photons by 
relativistic electrons, it is expected that the escaping 7-ray spectra from the binary 
systems are formed in IC e pair cascades. These cascades develop in the anisotropic 
radiation of the massive star since the injection place of primary particles lays 
outside the star. As a result of these cascading processes, a dip in the 7-ray spect rum 
is expected at energies determined by the temperature of the soft radiation field SSI. 
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4. General scenarios for the IC e^ 1 pair cascades 

Since the 7-ray signals observed from the massive binaries are modulated with their 
rotational periods, it is expected that 7-rays are mainly produced within the volume 
of the binary system where the medium is very inhomogeneous. Below we consider 
different effects which can influence the emerging 7-ray spectra. 



4.1. IC e ± pair cascades: linear, isotropized, driven by magnetic 
field 

As we mentioned above, the cascading processes very likely play an important role 
in the process of 7-ray production. However, they can occur in different scenarios. 
In the simplest case, the cascade develops through the IC and 7 — 7 absorption 
processes along the direction of the first generation of produced 7-rays (see Fig. 3 
on the left). We call this cascade as a linear cascade. Such approximation of the 
geometry of the cascade is valid provided that the Larmor radii of the secondary 
pairs, produced in the cascade, are much larger than the characteristic IC scattering 
length (the mean free path for the IC scattering) . This condition is usually not met 
within the considered binary systems due to relatively strong magnetic field. For 
detailed calculations of the IC pair cascade process in this scenario see Cerutti 
etal.Ha. 

In a more realistic IC pair cascade process, it is assumed that the secondary 
cascade pairs are completely isotropized in the place of their creation within the 
binary system. Again, the important criterion for the isotropization condition is the 
comparison of the Larmor radii of leptons with their IC scattering lengths (both 
dependent on particle energy and the magnetic field strength and the distance from 
the star). In this cascade scenario, the next generation of 7-rays is produced by 
pairs being isotropized in the place of their origin. These leptons produce 7-rays 
preferentially in the direction towards the source of the soft radiation field, i.e. the 
massive star (see Fig. 3 in the center). Some of them can also emerge from the 
binary system on the opposite side of the star with respect to the location of the 
source of primary 7-rays (or ele ctrons). Such type of the cascade has been at first 
investigated by Bednarek ^^ l 

In the third type of cascade (the most realistic one), secondary cascade pairs 
follow the direction of the local magnetic field lines at the place of their origin 
(see Fig. 3 on the right). In such a case, the secondary 7-rays can be produced in 
completely different directions than the direction of propagation of previous gener- 
ation of 7-rays. As a result, 7-rays, escaping from the binary system, can produce 
a characteristic pattern on the sky around the binary system which depend on the 
structure of the magnetic field around the star and the geometry of injection of 
primary electrons (i.e. the angle towards the center of the star, see Fig. 4). Such 
type of cascade has been at first studied by Sierpowska & Bednarek ^3. 
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Fig. 3. Inverse Compton pair cascades in the anisotropic radiation field of the massive star 
can occur in different geometrical scenarios. Left figure: The linear cascade in which secondary 
leptons and 7-rays move along a straight line (magnetic field is too low to disturb the directions 
of leptons); Middle figure: Secondary cascade leptons are locally completely isoptopizcd by the 
random component of the magnetic field in the stellar wind. Secondary 7-rays are produced pref- 
erentially in the direction towards the massive star. Right figure: Secondary leptons follow the 
direction of the local magnetic field lines (the most general case) producing the next generation of 
7-rays in a completely different direction than the direction of their parent 7-rays. 



4.2. Effects of e ± pair energy losses in the magnetic field 

The energy losses of primary electrons and secondary cascade e pairs on the syn- 
chro tron process can play important role on the development of the IC cascade 
(see | H|51 | 52 | 53 fy j n f ac ^ com p ar j n g the energy loss rates of electrons on the 
synchrotron and the IC processes, we can estimate the limit on the magnetic field 
strength below which the IC losses starts to dominate. When the ICS occurs in the 
Thomson regime, the limit on the magnetic field is simply given by the condition 
B < 40T| G, where T = 10 4 T4 K is the surface temperature of the companion 
star El F or typical surface temperatures of the stars in these binaries, of the or- 
der of 3 x 10 4 K, this critical value of the magnetic field is equal to ~ 400 G, i.e. 
not so restrictive in the main volume of the binary system. In the Klein-Nishina 
regime, the condition on the magnetic field strength becomes more restrictive H| 
Detailed calculations have shown that when the magnetic field in the cascade region 
is above ~ 1 G , then significant energy is transferred from leptons to synchrotron 
radiation E2E3] ^ g a resu [-j ;) jq e ± p an - cascade can start to be inefficient especially 
at the TeV energies when the IC cross section in the Klein-Nishina regime drops 
significantly. Note that the energy density of the magnetic field around the star 
drops in the interesting range of distances from its surface as Ub oc i? -4 , whereas 
the energy density of the stellar radiation drops according to Ub oc R~ 2 . There- 
fore, the region in which the synchrotron losses of electrons dominate over their IC 
losses should extend rather relatively close to the companion star. At the distance 
of 4 — 15 stellar radii from the companion star, the magnetic fields are likely to be 



November 28, 2011 1:20 WSPC/INSTRUCTION FILE Bednarek'l 



EXPECTED GAMMA-RAY EMISSION FROM X-RAY BINARIES 7 




Fig. 4. The distribution of 7-rays on the sky (darker regions mean more 7-rays) produced in the 
IC pair cascade in which the secondary leptons follow the local direction of the magnetic field 
lines (in the place of their origin), see for details Sierpowska & Bednarek (2005). The injection 
angles of primary 7-rays are equal to: 90° (upper panel), 120° (middle panel), and 150° (bottom 
panel) . 

of the order of ~1 G, in accordance with the calculations | H | 52 |53J 

4.3. Dependence on the shock localization (variable stellar wind) 

The 7-ray spectra emerging from the anisotropic IC pair cascade strongly depend 
on the shape of the compact object orbit. Therefore, the 7-ray spectra arriving to 
the observer should strongly depend on the injection place of electrons within the 
binary system which is in fact determined by the structure of the shock and also the 
orbital shape of the compact object. In Fig. 5, we show how these angles can change 
in the case of the shocks which bound around the massive star and the compact 
object. The shape of the shock depends on the pressure of the stellar wind which 
might change with time and the phase of the compact object (e.g. due to the effects 
of irradiation of the stellar surface or the anisotropic stellar winds). Therefore, we 
expect the the 7-ray light curves from specific binaries (especially those with variable 
winds, e.g. produced by Be stars) may change from one to another binary cycle. In 
fact, such time dependent 7-ray light curves have been observed in the case of LS I 
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Fig. 5. The angle between the direction to the observer and the site (a part of the shock region) 
where electrons are accelerated can change significantly (from the angle «i to 02), depending on 
the localization of the shock within the binary system. The location of the shock can shift as a 
result of the change of e.g. the stellar wind pressure. 

+61 303 which likely contain a Be type star. 

4.4. Effects of anisotropic stellar /pulsar winds 

As we mentioned above, winds of massive stars are expected to be strongly 
anisotropic. Howe ver also the wind produced by the pulsar is expected to be 
anisotropic I 54 | 55 [ Therefore, very complicated interaction patterns (shock struc- 
tures) between such anisotropic winds are expected (for some examples see Fig. 6). 
As a result, depending on relative geometrical orientations, the shock created in 
the wind coll isions can sometimes appear very close to the massive star or to the 
neutron star I 56 | 57[ j^- j s expected that such drastic changes may occur at quite 
nearby phases of the compact star on its orbit around the massive star. Thus, it is 
not very surprising that the 7-ray light curves from the binary systems can show 
very irregular features and high level of emission at unexpected phases. In specific 
phases the shock structure can sometimes appear very close to the star where the 
conditions for 7-ray production are the most favourable. 

4.5. Effects of clumpy stellar wind 

The winds produced by massive stars are expected to contain dense clumps im- 
mersed in a diluted medium filled with matter with different properties (magnetic 
field strength, density). The accelerated particles can diffuse into such two com- 
ponent medium suffering different propagation and radiation effects. As a result, 
the produced 7-ray spectrum may have a complicated structure due to the main 
energy losses on different radiation processes. Such model has been consid ered as 
a possible explanation of the emission features observed in LS I +61 303 In 
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Fig. 6. Possible geometries of the shock structures (shock I and shock II) which might appear 
within the binary system as a result of the interaction of axisymmetric pulsar wind with axisym- 
metric stellar wind which is composed of the equatorial dense and slow wind and poloidal fast and 
rare wind. 



the case of microquasars, the interaction of such a dense clumps with the particles 
accelerat ed in t he jet may also lead to the production of 7-rays through hadronic 
processes I 59 | 60 | 



4.6. Effects of relativistic boosting of radiation 

The relativistic pulsar wind, after slowing down in the downstream region of the 
shock, is still expected to move with substantial velocity, of the order of 0.3 — 0.5c. 
This is still enough to collimate at some level the 7-ray emission produced by 
relativistic particles in the wind 1 ^. In specific conditions, the relativistic flow after 
the shock is even expected to accelerate to very large Lorentz factors, of the order 
of ~ 100 In such a case, strong spikes of 7-ray emission should be observed in 
some phases of the orbit. Similar relativistic boosting effects, to those observed in 
active galaxies, ha s b een also consider ed w ith the application to microquasars by 
e.g., Romero et al. and Dubus et al. E3] g ucn an effects might be responsible for 
the modulation of the 7-ray signal with the period of the binary system. 

4.7. Double shock structure - acceleration of two populations of 
electrons 

In fact, the interaction of the winds from two stars turn to the production of a double 
shock structure separated by a contact discontinuity. As a result, the conditions 
for acceleration of particles on both shocks can differ significantly due to different 
fluids (characterised by different magnetic field strength, plasma velocities, etc ...). 
Therefore, it is expected that particles are accelerated to different maximum energies 
on both shock. As a result, two populations of electrons/hadrons are injected into the 
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binary system. Such a model can provide a nice explanation for the complicated (two 
component) 7-ray spectra observed from the binary system of two massive stars, Eta 
Carinae (see for details Bednarek & PabichGSl) j n the case of the binary systems 
containing pulsars, the shock at relativistic wind from th e p ulsar can accelerate 
electrons to maximum energies of the order of several TeV On the other hand, 
particles accelerated at the shock in the stellar wind can reach maximum energies 
of the order of a few tens GeV. Such a model provides natural explanation of 
the two component 7-ray spectra observed from LS 5039 and LS I +61 303 as a 
consequence of production of 7-rays by two populations of electrons with different 
maximum energies. Note that the kinetic powers of the winds from the massive stars 
in LS 5039 and LS I +61 303 seem to be comparable to the GeV 7-ray luminosity. 
However, the 7-rays produced in the IC e^ 1 pair cascades are highly anisotropic 
which result in enhanced 7-ray fluxes escaping in specific directions. 



5. Acceleration of electrons and/or hadrons ? 

It has been argued that leptons may have problems with reaching the TeV en- 
ergies within the binary systems due to the huge energy losses^. In such a case, 
hadrons can be responsible for GeV- TeV 7-rays observed from binary systems. Su ch 
a hadronic models have been considered in a few papers I 67 | 68 | 69 | 70 | 71 | 60 | 59 | 65 | j n 
fact, complicated 7-ray spectra observed from binary systems might be in agree- 
ment with the expectations of such models. However, the problem appears whether 
sufficient amount of energy can be transferred to radiation from hadrons due to 
their relatively inefficient cooling. The definitive confirmation of the importance 
of hadronic radiation processes should come from the observations of the neutrino 
signals from the 7-ray binaries. 



6. Conclusion 

Massive binary systems are one of the best defined sites for high energy processes 
and 7-ray production due to the well known geometry, soft radiation field, density of 
matter and structure of the magnetic field of the massive companion star. However, 
they are still quite complicated objects with complicated geometry. Due to the 
present lack of precise information, it is difficult to fix the basic conditions in these 
systems. Therefore, reliable predictions of the 7-ray emission features from 7-ray 
binary systems are difficult since many effects have to be taken into account such as: 
geometry of the particle interaction and the binary system in respect to the observer, 
non-steady medium (aspherical, variable, inhomogeneous winds), different radiation 
processes, different populations of accelerated particles. It is not surprising that 7- 
ray emission features from massive binaries show such variety of fine structures and 
unexpected behaviour. 



November 28, 2011 1:20 WSPC/INSTRUCTION FILE Bednarek'l 



EXPECTED GAMMA-RA Y EMISSION FROM X-RA Y BINARIES 11 



Acknowledgments 

This work has been supported by the Polish NCBiR grant: ERA-NET- 
ASPERA/01/10. 

References 

1. Bignami, G.F. et al. A&A 55, 155 (1977). 

2. Vestrand, W. T., Eichler, D. ApJ 261, 251 (1982). 

3. Maraschi, L. , Treves A. MNRAS 194, 1 (1981). 

4. Weekes, T.C. Sp.Sci.Rev. 59, 315 (1992). 

5. Protheroe, R.J., Stanev, T. ApJ 322, 838 (1987). 

6. Moskalenko, I. et al. MNRAS 260, 681 (1993). 

7. Paredes, J.M. et al. Science 288, 234 (2000). 

8. Mori, M. et al. ApJ 476 842 (1997). 

9. Thompson, D.J. et al. ApJS 101, 259 (1995). 

10. Vestrand, W. T. et al. ApJ 483, 49 (1997). 

11. Bednarek, W. A&A 278, 307 (1997). 

12. Bednarek, W. A&A 362, 646 (2000). 

13. Aharonian , F.A. et al. A&A 435, L17 (2005). 

14. Aharonian , F.A. et al. A&A 442, 1 (2005). 

15. Albert, J. et al. Science 312, 1771 (2006). 

16. Abdo, A. A. et al. ApJ 736, 11 (2011). 

17. Tarn, P.H.T. et al. ApJ 736, 10 (2011). 

18. Dhawan, et al. in Proc. VI Microquasar Workshop: Microquasars and Beyond. (Como, 
Italy, 2006) p. 52.1. 

19. Ribo, M. et al. A&A 481, 17 (2008). 

20. Albert, J. et al. ApJ 693, 303 (2009). 

21. Aharonian , F.A. et al. A&A 460, 743 (2006). 

22. Anderhub et al. ApJ 706, L26 (2009). 

23. Takahashi, T. et al. ApJ 697, 592 (2009). 

24. Abdo, A. A. et al. ApJ 706, 56 (2009). 

25. Abdo, A. A. et al. ApJ 701, 123 (2009). 

26. Acciari, V.A. et al. ApJ 738, 3 (2011). 

27. Jogler, T. et al. in Proc. 32nd ICRC (Beijing, China, 2011). 

28. Aharonian , F.A. et al. A&A 469, LI (2007). 

29. Acciari, V.A. et al. ApJ 698, L94 (2009). 

30. Corbet, R.H.D. et al. ATEL:3221 (2011). 

31. Tavani, M. et al. Nature 462, 620 (2009). 

32. Abdo, A. A. et al. Science 326, 1512 (2009). 

33. Aleksic, J. et al. ApJ 721, 843 (2010). 

34. Albert, J. et al. ApJ 665, L51 (2007). 

35. Sabatini, et al. ApJ 712, L10 (2010). 

36. Malzac, J. et al. A&A 492, 527 (2008). 

37. Del Monte, E. et al. A&A 520, 67 (2010). 

38. Acero, F. et al. A&A 508, 1135 (2009). 

39. Saito, T. et al. in Proc. 31st ICRC (Lodz, Poland) larXiv:0907.1017l f2009l. 

40. Aharonian , F.A. et al. A&A 439, 635 (2005). 

41. Hayashi, S. et al. APh 32, 112 (2009). 

42. Aleksic, J. et al. ApJ 735, L5 (2011). 

43. Tavani, M. et al. ApJ 698, 142 (2009). 



November 28, 2011 1:20 WSPC/INSTRUCTION FILE Bednarek/1 



12 Wlodek Bednarek 



44. Abdo, A. A. et al. ApJ 723, 649 (2010). 

45. Farmer, C. et al. A&A 526, 57 (2011). 

46. Walter, R. et al. in Proc. ICRC (Beijing, China, 2011). 

47. Usov, V.V., Melrose, D.B. ApJ 395, 575 (1992). 

48. Aharonian , FA. et al. J.Phys.Conf.Ser. 39, 408 (2006). 

49. Cerutti, B. et al. A&A 507, 1217 (2009). 

50. Sierpowska, A., Bednarek, W. MNRAS 356, 711 (2005). 

51. Bednarek, W. & Giovannelli, F. A&A 464, 437 (2007). 

52. Bosch-Ramon, V. et al. A&A 489, 21 (2008). 

53. Khangulyan, D. et al. MNRAS 383, 467 (2008). 

54. Bogovalov, S.V., Khangoulian, D.V. MNRAS 336, L53 (2002). 

55. Volpi, D. et al. A&A 485, 337 (2008). 

56. Sierpowska-Bartosik, A., Bednarek, W. MNRAS 385, 2279 (2008). 

57. Sierpowska-Bartosik, A., Torres, DF. ApJ 693, 1462 (2009). 

58. Zdziarski, A. A. et al. MNRAS 403, 1873 (2010). 

59. Owocki, S.P et al. ApJ 696, 690 (2009). 

60. Araudo, A.T. et al. A&A 503, 673 (2009). 

61. Dubus, G. et al. MNRAS 404, L55 (2010). 

62. Bogovalov, S.V. et al. MNRAS 387, 63 (2008). 

63. Romero, G.E. et al. A&A 393, L61 (2002). 

64. Dubus, G. et al. A&A 516, 18 (2010). 

65. Bednarek, W., Pabich, J. A&A 530, 49 (2011). 

66. Bednarek, W. MNRAS 418, L49 (2011). 

67. Romero, G.E. et al. A&A 410, 1 (2003). 

68. Romero, G.E. et al. ApJ 632, 1093 (2005). 

69. Kawachi, A. et al. ApJ 607, 949 (2004). 

70. Chernyakova, M. et al. MNRAS 372, 1585 (2006). 

71. Torres, D.F., Halzen, F. APh 27, 500 (2007). 



